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1 . INTRODUCTION 


Stall and its consequences are fundamentaUy important to the design and 
opera “i ofnisht vahlSlaa. A certain ■‘^^ree of unnteaaineas^ 

fin« nwr Art airfoil or other streamlined body at iiign angle 
attack, but the stall of a lifting surface undergoing unsteady motion is even 
more complex than static stall. Although much P^^^ress has bee^^^^ 
recent years, dynamic stall remains a major unsolved problem with a variety 
currenc^applicaUlons in aeronautics, hydrodynamics, and wind engineering. 

These lectures will summarize the main physical features of the phenomenon 
stall are discussed principally in physical terms. 


2. GENERAT. FEATURES 


Above a certain critical angle of attack, the flow around a slender lift- 
inn surface breaks down Into tlie phenomenon called stall. On an oscillating 
:i.ose”,.cldc„cc ic mere, .am. mcldl, u,.e 
delnved to incidences considerably in excess of the static stall angU . unce 
dynamic stall does occur, however, it is usually more severe and moic F^^sis- 
tLt than static stall. The attendant aerodynamic forces and moments exhibit 
large amounts of hysteresis with respect to the instatitaneous angle of attack 
a(t). especially if o oscillates around some mean value «o 
order of the static stall angle 0 ^,^.. 


2.1 Dymimic Stall Events 


Nnmeri'xis exper iments 
Have sliown that dynamic 
stall is cli iracteri/.i'd by 
llie shedding and p.issage 
over Llie upper surface ol 
I He lilting surface of a 
vei l ex-1 ike d i sluibanee . 
i'His viscous disturbance 
indmes a liiglily nonllnea ', 
Ituctu.Uinj; pri’ssma’ field 
(illust rated in I'ig. 1). 

If t He frequency, ampli- 
tude, and maximum inci- 
di'ucc ari’ sii 1 I i c i I’lil 1 y 
liij'.li, ti)c vortcx-.sbcddiuv, 
pHeuoru’Uon .s v.'c I 1 I'fg.a- 
n i .cii a nil c 1 oa t' 1 y d c t i lU'd . 



« ' *1 



iUuRtnIrod“?n passffl throuRh the di«tlnct stages 

t I' r‘Ta'.°"“ '^'■‘^ Table f ahnw 

Fig. 3 indicitos tlu> versus angle of attack, while 

1,1 dSagr*? ' ““ 


n the unseparatod region between points 1 and 2 in Fles 2 and 3 


llneir th^n Ti -i'P-*aavo...«i-ej.y une crenaa ot unsteady 

linear thin airfoil theory. Sometime after o exceeds a a ti^n li»Lr 
of reversed flow dPvoiT.„« r .«sa* a thin layer 


of reversed flow develops at the bottom of the boundary .'a?er. On the 
so-called trailing-edge stalling airfoils (e.g.. Fig. % this toL.o of 
5^“'^ Starts at the rear of the airfoil and moves forward to the 
ng-e ge region, whereas on leading-edge stalling profiles it develops 




TABLE 1.- THE DYNAMIC STAI.l, EVENTS OF FIGS. 2 AND 3 


Point 

Flow structure 

Forcea and momenta 

1 

Thin attnr.hed boundary layer 

Linear regime 

2 

Flow reversals within boundary 

Exceed static ^L^ax* 


layer 

linear regime 

3 

Vortex detaches and moves over 

Pitching mf .t diverges, vortex 


airfoil surface 

lift preset , 

4 

Vortex continues toward trailing 

Maximum lift and moment, followed 


edge 

by rapid decay 

5 

Secondary vortex 

Secondary peaks 

6 

RG£ittiicliiii6nt of flow fcom 
edge 

Readjust to linear regime 


quickly and very locally just dovmstream of the suction peak on the upper 
surface. In either case, a vortex then begins to evolve near the leading- 
edce region and spread rearward, as shown in the two center scones in Fig. 1, 

.,L a speed somewhat less than 1/21U. The associated distortion in the 
pressure distribution causes the quarter-chord pitching moment to diverge 
from its previous trend, at point 3 in Figs. 2 and 3, to large negative 
values. The drag also begins to rise dramatically. However, unlike static 
stall, the lift usually continues to increase mouotonically until the vortex 
is well past midchord. 

As the vortex nears the trailing edge, lift, moment, and drag reach their 
largest values (point 4). although usually not exactly simultaneously, and 
then drop dramatically. Secondary and sometimes tertiary vortices pioduce 
additional fluctuation in the airloads, point 5. but at greatly reduced ^<^vels. 
If a is decreasing in the meantime, the flow will begin to reattach, point 6, 
at some angle of attack much less than The reattacliment point moves 

rearward at a speed well below so that several chord lengths of travel 

are required before the flow completely returns to the approximately linear 

domain. 

The dynamic airloads are, of course, manifestations of the unsteady pres- 
sure distributions on the airfoil. Therefore, it is Instructive to look at 

as a function of space .and time. Figure 4 shows the pressure distribution 
ai discrete times, obtained from lb miniature pressure transducers distributed 
over the upper surface of the airfoil. Figure 3 s'lows the continuous time 
histories of each of these 16 tr.ansducers . The abrupt collapse of suction 
around the leading edge and passage of tlie primary and secondary vortices i c 
partlcvilar ly evident in these figures. 


2.2 Hvsteresis .ind AertHlyiuunic Damping 


bach of the aforementioned dvv.amlc stall events takes a finite .umnmt ot 
time t'o develop, scaling approximately as U^.At/c. However, once they are 
Initiated, thel Lend to he relatively Independent of the airfoil motion. In 
addition to the obvious implications for mathematical mode ing of the phenom- 
U,l» h,.s cwo co„s,-q.u.n.n.». The irst I» « -.J 

asymmetry of the airloads with respect to the motion o the ^"’dy ih is P o ^ 
duces (he hysteresis discussed in Section 2.1 and read I y app.uent in Ig. 

I his vei'iesents an Imporlanl eoilrast to the quas i -sieaily limit, fot which the 
How field adiusis immediately to e.ich change in incidence. 


In the steady 


- t 


.60 



Fiij. 4. Pvetisuve (Hsti^i-hutions at 
times dwhu^ dyuamie stall 
fi'}' the eonditions of Fig. 2. 


rc "°®® PRIMARY SECONDARY 

VORTEX VORTEX 

t 1 ' i 

-W2 0 nl2 ’ 

Fig. 5. Tima historn.es of pressures 
at var-Lous poiaits on the air foil of 
Fig. 2. 
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case the aerodynamic forces are approximately unique functions of a. whether 
a is increasing or decreasing; therefore, angle of attack becomes a primary 
parameter to the aerodynamicist . It remains an important one in unsteady 
separated flow, but it is clearly not so unique. 

The second and closely related point concerns what the aerodynamicist 
considers to be aerodynamic damping or. more properly, the net aerodynamic 
vVorfoer cycle of oscillation, and its relationship L- flutter For a pitch- 
ing airfoil, the instantaneous work done on the fluid by the body due to its 
motion is dW - -Man, where M is the pitching moment about the axis of rota 
tlon and is positive for nose-dowiT airloads. Although the product of i-Mdu 
is normally negative, during some phases of dynamic stall it can become 
positive, so the fluid is doing work on the body instead of vice 

The so-called pitch-damping parameter, given by f. -|c„du/4aiS is a 
measure of the net work done over a complete cycle ol oscillation. If t is 
neaatlve the airfoil extracts energy from the airstream. and the pitch osc 
laaottr^iVi tend to increase in .amplitude unless restrained. This, o course 
is the condition for flutter, and the hysteresis ol dynamic stall pemitb it 
to occur in a single degree of freedom of oscillatory body motion. Normally, 
in unseparated I low.s f hitler only occurs when tlie body mot. on iiKludcb 
multiple degrees of freedom, e.g. , combined bending and torsion of an aircraft 

wing. 

Stall flutter, arising from this negative pitch damping, tends to occur 
when the airfoil is oscillating in and out of stall. As Indicated in Hg. 6 
from Ref. 10. the damping is given by the area inside the 0^,-0 trace. i 

Lh's example, it is positive throughout the cycle when . i-t-., w‘hu 

1.0. . wlu-.. 0,o airfoil Vo., a.no at.nllod 

Lh-oughouL the cvcle. However, for = 14.9". the airfoil is sl^illcd part 

Iho .1..0 a„.l :,ot at oil, or tlmoa. and tho area,. Inaldo t l,o clockw.so 
(• •. 0) and countercloekwise loops (r, > 0) are approximately the same. Ihis 

indicates neutral stahlUtv. Had the mean angle been slightly less, the ret 

hn, non U1 „.ol,al,lv l.avo l.oon noRUtlvo ,n„l tUo oncll lat Ion nonUI Iho.oroio 
leiv.- been unsl able if uni e.str.i i ned . 


M„‘i0.40 a Off tin cot 
STATIC DATA 


POSITIVE DAMPING 


P NEGATIVE DAMPING 



Stnll-lnducQd negative damping 
can occur In other types of motion, 
ouch aa plunging OBcllIatlona 
(Ref. 10 ). In this case, the aerody- 
namic work is dW - Ldh, where h is 
the displacement of the airfoil per- 
pendicular to the oncoming flow. As 
before, unstalled motion is normally 
stable, and the tendency toward 
instability is greatest when the air- 
foil oscillates in and out of stall. 


3 . STALL REGIMES 


The flow field around an oscil- o 

lating airfoil in subsonic flow can 
be characterized by the degree or , 

extent of flow separation. For a 
given airfoil, the primary parameter 
that determines the degree of separa- o ^ 

tion is the maximum angle of attack c 
(”max “ ‘^o + for sinusoidal 7.3 • 

oscillations). An important aspect ~ 

of dynamic stall is the large ampli- 0 'o ''<n 14.9 ^ 

tudcs of the motion which produce , 1 L I | 

the large maximum angles. This con- *0 10 20 

trasts with the hierarchy of viscous incidence, a, deg 

effects on oscillating airfoils at ri* . , - 

transonic speeds and low angles of cjfect of mean < 

attack, where the scale of the inter- moment ooefftoients 

action is governed primarily by the aei'odyyiamio damping. 

strength and motion of the shock 

wave. So far, prediction methods for this class of problems have not 
successful for the low-speed, high-angle problems and vice versa. 


Fig, 6. The effect of mea>i c 
Tift and moment coefficients 
aemduriamio damping. 


The Importance of fs illustrated in Fig. 7 (from Ref. 8), which 

portrays four important regimes of viscous-lnviscld interaction for oscillat- 
ing airfoils. For the left-hand part of the figure, i.e., = 13 % there 

is almost no separation throughout the cycle, although unsteady boundary-laye 
displacement thickness effects are not completely negligible, t^hen is 

increased to U”, the limited separation chat occurs during a small fraction 
of the cycle distorts the hysteresis loops of the unsteady pressures and air- 
loads. From a practical standpoint, the effect on Cj^j is particularly impor 
taut. Tills stall-onset condition represents the limiting case of the maximum 
lift that can be obtained with no significant penalty in pitching moment or 
drag. 

A slight additional increase In Omiix produces a major increase 

In the extent, severity, and duration of the separation phenomenon, for the 
conditions shown In Fig. 7 . This type of viscous-lnviscld interaction pro- 
duces what is called "light dynamic stall" (Refs. 8-9). Further increases in 
‘huax to the deep dynamic stall regime, with ;i large viscous zone over 

the entire upper surface of the airfoil during half or more of the cycle. 
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Fig, 7. Dynamia stall regimes. NACA 0022 airfoil, a ^ + 10° sin wt, 

k = 0.10s solid lines denote increasing a, dashed lines decreasing a. 


3.1 Light Stall 

This category of dynamic stall shares some of the general features of 
classical static stall, such as loss of lift and significant increases in drag 
and nose-down pitching moment compared to the theoretical inviscid values , 
when « exceeds a certain value. In addition, the unsteady stall behavior is 
characterized by large phase lags and hystereses in the separation and 
reattachraent of the viscous flow and, consequently, in the airloads. Also, 
the tendency toward negative aerodynamic damping, as discussed above, is 
strongest in this regime. 

Another distinguishing feature of the light dynamic stall is the scale of 
the interaction. The vertical extent of the viscous zone tends to remain on 
the order of the airfoil thickness, generally less than for static stall. 
Consequently, this class of oscillating airfoil problems should be within the 
scope of zonal methods or thin- layer Navier— Stokes calculations with rela- 
tively straightforward turbalence modeling. 

The quantit.ativc behavir>r of light stall Is knowii to be especially sensi- 
tive to airfoil geometry, reduced frequency, maximum Incidence, and Mach 
number; also, three-dimensional effects and the type of motion are probably 
important. T!ie qualitative behavior is closely related to the boundary- 
layer-separation characteristics, e.g.i leading-edge versus trailing-cdge 
separation, and to the changes in this separation behavior with u^ax* 

The effects of these parameters are discussed in Section 4. 


3.2 Deep Stall 

The vortex-shedding phenomenon discussed previously Is the predominant 
charactei isLlc of this dynamic stall regime. The passage of the vortex over 
the upper surfai’c o\' the .lir'^oil produces values of Cjj, and Cp far In 

exces.s of their static counterparts when i is increasing, .and large amounts 
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of hysteresis occur during the rest of 
Che cycle. Ihe scale of the viscous 
interaction zone is also large; the 
thicknoss of the viscous layer is on 
the order of the airfoil chord during 
the vortex-shedding process. 

Figures 7 and 8 (from Ref, 8) 

Ulus crate some of the qualitative and 

quantitative differences in light and 

deep dynamic stall. Moment stall, 
denoted by M, occurs rather abruptly 
for all three values of Uq, but the 
deep stall drop in lift after C T.may i 
denoted by L, is not evident in the 
light stall case. The large negative 
values of are due to the vortex. 

The qualitative features of deep 
stall are less sensitive to Che details 
of the airfoil motion, airfoil geome- 
try, Reynolds number, and Mach number, 
provided strong shock waves do not 
develop (Ref. 8). The quantitative 
airloads depend primarily on the time 
history of the angle of attack for 
the portion of the cycle when a 
exceeds the static stall angle Ogg. 
This feature and other details of the f 
the following discussions. 


NACA 0012 Mf, o 04 u ^ + 10* tin fut 


MQMTSTAU 



M 




1 

1 1 


DEEP STALL 



Fiff. S, Effect of mean angle on the 
time hiatovy of pitching moment 
coefficient for the conditiona of 
Fig. ?. 

' behavior will become evident in 


4. EFFECTS OF VARIOUS PARAMETERS 


One of the reasons that dynamic stall is more difficult to analyze and 
predict than static stall is its dependence on a much larger number of param- 
eters. Table 2 gives a general indication of some of the more important ones 
and their effects. 


4. 1 Airfoil Geometry 

Especially in the light stall regime, the leading-edge geometry of an 
airfoil is a principal factor in determining the boundary- layer separation 
characteristics. Airfoils with moderately sharp leading edges tend to develop 
severe adverse pressure gradients in the first few percent of chord, leading 
to abrupt boundary-layer separation there that spreads rapidly downstream. 

This is referred to as "leading-edge stall," and it produces a relatively 
concentrated vortex and abrupt changes in Cj.j, and Cp during the develop- 

ment of the dynamic stall events. This contrasts with traillng-edge stall, 
which is more common on airfoils with rel.«tlvely blunt noses or large amounts 
of leading-edge camber. In this case, the boundary-layer separation progresses 
forwartl from the trailing edge, the onset of stall tends to occur more grad- 
ually, and unsteady effects tend to suppress the separation more than on 
leading-edge stalling airfoils. Also, the more gradual the traillng-edge 
separation, the less likely negative aerodynamic damping will occur. 

I'lg.ures 9 and 10 sliow some of the airfoil sections that have been 
studied by the author and his col leagiu^M (Refs. b-8). Results In the light 
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TABLR 2 . - im’ORTANCE OF DYNW'IIC STALI. 
PARAMETERS 


Stall parameter 

Effect 

Airfoil shape 

Large in some cases 

Mach number 

Small below - 0,2, 

Large above M^^ ^ 0,2 

Reynolds number 

Small (?) at low Mach number, 
Unknoitm at high Mach number 

Reduced frequency 

Large 

Mean angle, amplitude 

Large 

Type of motion 

Virtually unknown 

Three-dimensional 

effects 

Virtually unknown 

Tunnel effects 

Virtually unknown 



eiudied hi Rcfv. (> and ?, 


C. 


NACA0012 


AMtSOI 


C__ 




HUGHES HH 02 


VFRTOL VR 7 


WO HIM ANN FX098 


NLR 1 




SIKORSKY SC 1095 




NLR7301 

Fig. 10. AivfoitL' studied hi Hef, 8. 


stall and stall onset roRimes for the profiles in Fig. 10 are shoim in 
FIr, 11, which illustrates how widely the aerodynamic coefficients can vary 
for one set of unsteady conditions. 

The type of boundary-layer separation is listed for each section in 
FIr. 11. In this example, tile maximum incidence, 15", was more than 1" greater 
than the static stall angle of all the airfoils, except the NI.R-7301 . Conse- 
tiueully, the oilier seven sections exhibited some dynamic overshoot of f-Rmax* 
However, this came at the expense of large pitching moments on the. NACA 0012, 
Wortmann FX-0‘)8, Sikorsky SC-1095, Hughes HH-02, and N1,R-1 airfoils, and 
negative pitch damping on the 0012 and SC-1095. 


Tlie differences between these airfoils are diminished in deep stall, 
where tlie vortex shedding becomes fully developed. Figure 12 (from Ref. 8) 
.shows the lift, moment, and drag results for four of the eight sections in 
I'lg.s. 10 and 11. The hysteresis loops for the Ames-01, VR-7, and N1,R-1 air- 
foils are remarkably similar, diffeilng princinally in the angle for the 
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onaat of stall, In the magnitude of the peak forces and momenta, and In the 
strength of the secondary vortex at or near The same could be said 

for the other four nirfollB (not shown). 


4.2 Reduced Frequency 

An important parameter affecting dynamic stall is obviously the reduced 
frequency of the oscillation, but its influence depends on the stall regime 
and the type of boundary- layer separation. For example. Figs. 13 and 14 

opposite trends in the light stall regime for leading-edge 
and trailing edge stalling airfoils tested under identical conditions. On the 
other hand. Fig. 15 (from Ref. 2) is more representative of the trends in deep 
dynamic stall. In this case, the vortex-shedding phenomenon developed for 



Fig. 23. The effect of reduced frequencij oil the VR-7 airfoil at M = 0 30 
end <x ~ 20" t sin ...t. 



i’-lO 


♦I a 16' + 1(T 


"" STATIC □ POSITIVE DAMPING NEGATIVE DAMPING 



Fig. lb. The effect of reduced frequency in the deep stall regime. NACA 0012^ 
= O.IQ. 


k > 0,05 and the strength of the vortex appeared to become independent of 
reduced frequency for k > 0.15. 


4.3 Amplitude and Mean Angle 

For periodic motion, and cannot be completely separated since 
“max “ “o very important in determining the amount of separation. 

This was discussed previously in connection with Figs. 6 to 8; but in those 
examples, remained constant while varied. On the other hand. 

Fig. 16 shows the effect of varying a, with ct„ and the pitch rate parameter 
‘Stiax constant (o^^jj^c/l’,., = 2tjk).‘ Kor the deep stall cases of = 10° 

and 14'*, the vortex-shedding phenomenon begins while a is still increasing. 
The subsequent flow-field development is very similar for these two cases, 
although the different stall events occur at different phases in the cycle. 
However, when “ 6“, dynamic stall does not begin until after <-imax> 
occurs as u is decreasing. As a result, the vortex effect is less pro- 
nounced, and the results are more representative of light stall. This is not 
unexpected and, in the limit of very small amplitude, the vortex would essen- 
tially disappear. 

An interesting comparison can be made if combinations of u^, Uj , and k 
arc selected to match the a(t) histories over the portion of the cycle where 
.1 exceeds the static stall angle. One such example is sho\m ir. Fig. 17. In 
these two data sots, the pitch rates, a, are very nearly the same at 
a 1 7 , v/hlch is about where moment stall begins. The lift and moment data 
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STATIC DATA LJ POSITIVE DAMPING 


NEGATIVE DAMPING 


+ 0" iln (Jt 


k ° 0.24 


/( ^ STATIC 


rt* 16' ■* 10' ilnwt 


a " 16' + 14 Iln wl 








20 30 


INCIDENCE, a . deg 


I 1 1 J I 

10 20 3 


Fig. 16. The effect of a^lityide^ ^mas: = 0.06. 


CJ 10 ♦ lO’jmwt; k 0.10 

rt = 15 -i- 5 iin wt; k =0 15 


Fig. 17. Funamc tUPioaite fov ~ " 0,002.. 


are almost identical on the upstroke and throughout most of the stall events, 
whether viewed versus Ucot/c or versus ci. This would not be the case if the 
i(t) histories were not so closely matched. 

A close exiimination of numerous cases fmin Ref. 8 and elsewhere indicates 
tlie general rule that the better the match between the u(t) histories between 



“static stall* S>ax* “reattachraant » better the match 

dynamic coefficients, at least for pitching motions In deep stall. However, 
the extent to which this conclusion can be extended to light stall or to other 
types of motion remains to be established. 


4.-i Mach Number 

Even at relatively low free-stream Mach number, supersonic flow can 
develop near the leading edge of an airfoil at high incidence. This la shown 
schematically in Pig. 18, in comparison with the more conventional transonic 
airfoil problem. As M„ increases, eventually the supersonic zone would be 
expected to terminate in a shock wave of sufficient strength to cause boundary- 
layer separation and leading-edge stall. The impact of such shock- induced 
separation on dynamic stall in general, and on the vortex-shedding phenomenon 
in particular, has not yet been assessed to any significant extent. 



Figure 19 shows some effects of increasing Mach number on dynamic stall. 
The data are taken from Liiva et al. (Ref. 10); here progressively smaller 
mean angles were selected because of the decrease in the static stall angle 
with increasing Mach number. The similarity of the static stall characteris- 
tics at M =0.2 and 0.4 suggests that transonic shock-wave formation docs 
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not ploy a tola In either cnee, but ‘he etetlc « ^ °y. 0.6 

evidence o£ahock-l^^eeper« aevelopmont o£ 

BUggeat that the formation or an voot-tttea of the phenomenon remain. 


4.5 Other Types of Motion 

Mn«i- nf the available Information concerning dynamic Gtall has 

obtairef ofeSfoiu oaclUetinB ^n- 

quarter-chord. However. with moat o£ the resulta outlined 

periodic motion shown qualitative agreement wx 

above . 

Plunging oaclllatlons. or vertical translational 
equivalent fluctuating angle Laume that the resul- 

mean incidence Oq. Current engine , ^ those due to pitching oscilla- 

taut flow fields and airloads arc equivalent to S Fig! 20 (from 

tlons This appears reasonable for deep stall, a ' , _ 14 v,^ c^nii 

Rff U)! but significant differences have been reported for light stall 

_ m ^ t K 


(Refs. 12, 13). 


“off ' ^ lincjl, k - 0.15 



30 


F?V< Ccini\vvGi-»i of piioh tvhl 

phOUU' I’t'.fJf//.' i tfl't’J' i'tO: /• 


The studies of Maresca et al. 

(Ref. 11) also established the exis- 
tence of the vortex-shedding phenom- 
enon and strong hysteresis effects 
due to translational oscillations in 
the streamwise direction and at 
various oblique angles rel, ;ive to 
the oncoming flow. In the 

longitudinal-oscillation case, for 

example, unsteady effects were found 
to be strong enough to unstall the 
airfoil during part of the cycle, 
even when the incidence was consid- 
erably above the static stall angle. 

4.6 Three-Dimensional Effects 

Despite the fact that almost all 
practical devices develop three- 
dimensional flow, common engineering 
practice is to use so-called strip 
theory and the two-dimensional stall 
characteristics outlined above. 

This Is partly due to simplicity and 
partly due to the acute lack of three- 
dimensional information. 


.. Cl HiliJi-e et al. (Ref. 14 ) examined the effect of 

.,„..vp:ru “urritSiug 

Uy cc.mp...«nt norm,.! to kt.-.ll nml to ro.tuco tl.o 

-an hop, ( ns! ,t a.^o n.0 

iiiagiiiuide of the hysteresis loops. 
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PITCH AXIS 

Fig, 21. Hketifh of osoittaPing iwcpt iHtig expovimcnt 
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Fitmlly, the moaauremonta of McCroakey ai^d Flahcr (Ref, 15) cm a model 
hollcoptor rotor may be raontioned. At a apanwiue station r, the rotor blade 
olomont experiences a periodic velocity component Vj^ « ilv + V„ sin fit normal 
to the loading edge and a apanwlae component cos nt due to combined rota- 
tional and translational motion of the rotor. In addition, the angle of 
attack changes periodically due to a complicated combination of pitching and 
plunging motion. 

Figure 23 shows two cases from Ref. 15, where Y - rtt. For the deep stall 
case, a detailed analysis of the pressure data Indlcoted that the stall events 
and the airloads were essentially the same as those on an airfoil oscillating 
in pitch for 90 < H* < 270", provided the large-amplitude oscillations in a 
were approximately matched. The vortex-shedding phenomenon appeared at 
T s 210", followed by large transient overshoots in and that excited 
stall flutter at the torsional natural frequency of the blades. 


V^/ilR ■> 0.3B t/R - 0.76 


V.4 



I, . rf«g 

Fi< 7 , S3, Blade cloKent or.vh'oy^nent and airloads oti a model votov. 


However, these and other model rotor data indicate that three-dimensional 
effects produce significantly larger lift for 270" < T < 360, after the 
vortex passes off the blade and into the wake. At the moment, therefore, the 
use of two-dimensional airfoil data to predict the airloads on devices such as 
helicopter rotors and wind turbines should be viewed with a certain ivnount of 
caution. 


5. METHODS OF CALCULATION 


The physical aspects of the phenomenon of dynamic stall have been out- 
lined in previous sections. We consider now a number of special techniques 
that have been suggested for predicting the engineering quantities of interest, 
such as the instantaneous v'alues of lift, drag, and pitching moment. It should 
be emphasized that all of these methods are still being developed, refined, and 
improved, even in two spatial dimensions. Furthermore, all of the techniques 
invoke some sort of simplifying assumptions and approximations and are tailored 
to the specific features of some particular stall regime or to a relatively 
narrow range of parametora. 
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5,1 Discrete Potential Vortex Approach 


A nromlainR analytical approach to the deep dynamic stall roRlrao in 
I ^ ^ a'kir, finu tnkufl its cue from the discrete vortex model that has been 

sKoteted la Fig. 24 (a). la the alaplar ca„e 
of a cirLlar cylinder in the low subcritlcal Reynolds number regime, alternate 
or staggered rows of potential vortices produce the essential features of ^he 
?Lcef ^the cyllnde?. but the vortex spacing and frequency have to be speci- 
fied empirically. 
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i\). :'i. The dioevcie iv i-.'.’.v 


The vortex-shedding phenomenon on a thin flat 
1 series of emitted vortices by Ham (Ref. lo) , ns shown on the right in 
Fig. 24. Each vortex moves under the Influence ol all the othcis . t 

result is a tendency for the Individual filaments to a structure 

that resembles the experimentally observed features of dynamic stall. 


As In the case of the cylinder, assumptions must be made regardlttg the 
geometry and strength of the '^'J^tex emlssl^ .,'"lgneS\o eac^vor^ and in 

r:::ha!:ir^:;r^i:Li?nfX^bL^r^ the boundary-layor- 

♦ ^ LM ist ics ou the body. Ham started the process at ai 

7i;' r “uv ss;» .y nc .aid ,.d)Ltod the atrength of the vot tox o»la- 

' o iLftratogaation golota at the too odgoo of Iho. ploto. ’^la aaouott. 

to tapooiog a apocltl Kutto condition at tho loading odgo an well a» at tho 

trailing edge. 
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luimei ical potent ial-f low lechnUiue of Mesing (Rel . IHK 


Further refinement:) and extensions of the diserete vortex approach have 
1 I Iv Oro et al (Kef. l'‘l. Ihilike provloii.s analyses, the all loll 

:::;;;;:,,::^i!:vort i- 




are choaou to eatlafy the no-allp condition on tl.fi surface. Tho surface 
vortex layer ia divided into n small number of finite-length elementa, with 
one vortex filament just off the surface for each element, as showii schemati- 
cally in Fig, 25. At time step N, tho vortex induces a surface 

velocity AqN.k » [<i.N,k+i _ ^N,k]/[sk+i _ gk] that cancels the velocity that 
would otherwise exist for irrotational flow around the body. Each vortex 
filament is then allowed to convect to a new position at the next time step 
under the influence of the flow field due to the body and all of the 
vortices. As time progresses and new vortices are created near the body, 
"separation" and the clustering of vortices into something resembling the 
vortex-shedding phenomenon observed in experiments evolves naturally. 

Figure 26 shows a typical pattern of vortices and uhe resultant streamlines. 



NO SUP BOUNDARY CONDITION: 

(u - i| • - (v - v) ■ ■■■ 0 

SATISFIED BY CONDITION IMPOSED ON 
AT."'. 





Af'pi u\il i'/* i )u' 


■/ to 


.,v: 


oti’uj oii'fot/. 


In principle, the only free par;imeters in this method are the size of the 
lime steps and the number of elements used to represent the boundary layer 
around tlie body. However, since the number of generated vortices continually 
Ineroases, computational limitations dictated two further approximations in the 
calculations that have been performed to date (Ref. 19). The first was to 
allow clusters of vortices far away from the body to be replaced by a single 
eqtiivalent vortex. Second, the number of elements represont ing the boundary 
layer was severely limited after the first time step. The results so obtained 
appear qualitatively correct for dynamic stall at low Reynolds numbers, but 
the method needs further development to become a practical engineering tool. 


5.2 Zonal Methods 


Kspcclal ly 
relatively thin, 
can be modeled s 


In the light stall regliiu’, the viscous region teiuls ti’ remain 
For flow fields of this type (Fig, 27), the various regions 
eparately and coupled in some appropriate fashion. 
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Thln-bouniiary- 1 OVLT approach- The. classical Prandtl boundary- 
layer equations, Incorporating some sort of turbulence modeling and modified 
ti> Include unsteady terms, have been used by several authors (see Ref. 1) to 
investigate some of the cliaracterist Ics of dynamic stall. These studies have 
shed light on the delay in the onset of dynamic stall and on tlu’ differences 
between steady and vmsteady separation; but lacking any coupling between the 
viscous and inviscid regions, this approach has very limited quantitative 
usefulness. Also, it gives few clues to the mechanisms responsible for the 
vortex-shedding phenomenon, which is such a prominent feature of dynamic 
stal 1 . 

5. .1.2 Coup led vis cmts-in vlscid inter act ions- The analysis of Crlmi and 
Reeves” (Ri^f..’ 01 retained boundary - 1 aver concepts wherever possible, while 
allowing, a stfi'iig couvillng to exist lu’tween tlii' viscous flow in tlu; separtited 



zone and the aurroundlng Inviscld flow. These authors modeled the Inviscld 

J unsteady thln~alrfoll theory, representlnR the slr- 

Hff'’ distributed source and vortex singularities, 

s of the unsteady boundary-layer equations were 

l erfomed foi the attached . low, using an eddy-viscosity model for the turbu- 
lent flow. In the absence of boundary-layer separation, no interaction 
between the viscous and invlscid flow was considered. The leading-edge sepa- 
ration bubble formation was analyzed to determine whether reattachment would 
occur or whether the bubble would "burst," and this criterion was used to 
distinguish between leading-edge and tr ailing-edge stall. 

A number of approximations and assumptions were made in developing the 
analysis, but most of the essential flow elements seem to have been included. 
Although the application to tralling-edge stall was not attempted, the method 
was applied to a variety of unsteady airfoil and helicopter problems, with 
varying degrees of success. One of the main faults was that the basic pre- 
diction of whether the leading-edge bubble would burst did not seem to depend 
on Reynolds number and leading-edge geometry in the proper manner. Attempts 
by the present author to apply the method to the prediction of the static and 
dynamic stall characteristics of .several helicopter-type airfoils which are 
thought to stall by the bubble bursting mechanism did not correlate well with 
experiments. Furthermore, the viscous part of the analysis does not account 
for the feature of a thin layer of reversed flow near the wall before and dur- 
ing the beginning of tralling-edge stall, nor for highly organized vortex- 
shedding, as discussed earlier. 

More recently, Rao et al. (Ref. 21) extended a steady viscous-tnviscid 
interaction method (Ref. 22) to d>Tiamic stall by introducing an empirical 
delay parameter. At, Into an effective angle of attack, 

'^eff = sln(,.‘t - .oAt). The parameter Ai is considered to be a func- 

tion of v^j and k tor a given airfoil and is constant througliout the cycle. 
The potential flow at each value of is calculated by panel methods, 

with a distribution of sources and vortices representing the airfoil and 
vortex sheets representing free-shear layers. The free-shear layers approxi- 
mately coincide with the edges of the viscous lavers in Fig. 27 and arc sllp- 
s u-face streamlines. They define an effective wake, which is assumed to\e 
liiotatioiyil and to have a constant total pressure less than that of the free 
struiun. The point where the upper-surface vortex sheet leaves the airfoil is 
determined by an integral boundary-layer calculation of the separation point 
(point S in Fig. 27). The resultant configuration is .sketched in Fig. 28. 
Iterations are performed at each value of a^ff: until the wake shape and’ 

sepaiatlon points converge; actually, the convergence criterion is tliat Cr 
change le.ss t!un IT.:. *-• 
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Although realistic hyatoroais loops of and Cjj voi-rub a have been 
cnleulotod, the results published so far (Ref. 21) cannot bo considered pre- 
dictive because they depend heavily on the arbitrary choice of At. And as in 
the analysis of Crimi and Reeves (Ref. 20), neither the special features of 
unsteady reversod flow nor strong vortex shedding have been included. However, 
with additionol modifications and refinements, these zonal methods could offer 
a reasonable balance of rigor, accuracy, and computational efficiency, at 
least for light dynamic stall. 


5.3 Navier-Stokes Calculations 


The limitations and approximations of zonal methods aiid the questions of 
unsteady separation models can, in principle, be avoided by solving the full 
Navier-Stokes equations. However, turbulence must be modeled and a number of 
computational difficulties must be overcome if this approach is to be realistic 
for most practical applications. 

Several recently developed methods for laminar flows, discussed in Refs. 1 
to 4, have shed considerable light on the physical mechanisms of deep dynamic 
stall at low Reynolds numbers and have laid the computational foundations for 
hlgh-Reyno Ids -lumber solutions of the so-called Reynolds-averaged, Navier- 
Stokes equations. For incompressible flow in two dimensions, they may be 
expressed as: 

V«q = n (O 

.’TTj u 

+ q*Vq = - — Vp + vV'q - — r— — (2) 

^ ^ P 

where q represents the veiociLy vector (u,v) time -aver aged over an interval 
sliort with respect to the time scale of t he mo tion, but long with respect to 
the time scale of the turbulence, tind -puj^ui is the Reynolds stress due to 
the tvindom turbulent fluctuations. The Reynolds stress term vanishes in 
l;minar flow and is to be modeled (empirically or otherwise) in the turbulent 
case. The particular representation of the detailed physics of turbulence is 
crucial in many steady-flow problems that are dominated by viscous effects. 
However, the importance of turbulence modeling In dynamic stall has yet to be 
established . 


Numerical codes using various turbulence models that have been developed 
for steady flow arc currently under development at several research centers. 

A few very preliminary results have been published (Refs. 23, 24) and more 
will become avail.ible in the near future. However, quite apart from the ques- 
tion of the large computational times and storages required, it will be some 
time before these ci'des are .adequately tested and verified. Nevertheless, the 
careful ;md selective analysis of new numerical results for both light and 
deep dynamic stall will be especially beneficial in guiding the development of 
more approximate engineering analyses and prediction methods basod on empirical 
correlat ions of wliul-tuimel data. 



S,4 I-juplrlcal Cone I at I on Techniques 

The helicopter industry has developed several engineering prediction 
eehnlques based on empirie.il eorielations i>t wind-tunnel data fi>r estimating 
he unsteadv airloads on oscillating alrtoils. These methods seek to corro- 
.ite force and moment data obtained from rolativelv simple wind-tuimol tests 
n fonrnil at ions that slu'W the effects o1 the niimei oiis relevant p.iramotevs 


-J1 


such ns airfoU shape , Mach number, amplitude and frequency of osciUntions, 
mean angle, and type of motion. 

Common to all the available literature relevant to dynamic stall Is the 
observation that unsteady effects increase with increasing pitch rate, that 
la, rate of change of airfoil incidence. It la also evident that the dynamic 
stall events require finite times to develop. Therefore, sonie form of the 
dimensionless parameters ac/U„ and U„At/c appears in all of the empirical 
methods. Another common aspect is that the empirical correlations are used as 
corrections to steady airfoil data, so that most geometrical, Reynolds niunber, 
and Mach number effects are only accounted for insofar as they determine the 
static section characteristics. 


The highlights of several methods currently in use by the helicopter 
Industry are outlined below. More detailed accounts of each can be found in 
the references cited and in Refs. 2 to 4, 


5.4.1 Boeing-Vertol gamma function method (Refs. 25, 26)- The onset of 
dynamic stall is assumed to occ ur at njjg = agg + Aup, where agg is the 
static stall angle and Aop = /U„. The quantity y» which is the essen- 

tial empirical function, depends on airfoil geometry and Mach number and is 
different for lift and moment stall. The gamma functions were generated from 
a large amount of data generated in a transonic wind-tunnel test of various 
airfoils oscillating sinusoidally in pitch. The force and moment coefficients 
are constructed from static data using an equivalent angle of attack that 
accounts fo r unstea dy potential-flow effects, ci , and a reference angle, 
oir = « ± >* follows: 


C = [a /a )C (a ) 
L eq L' r 


S “ S^«r) 


/ based on 


lift 


(3) 


Si “ 

The location of the aerodynamic center of pressure, is specified empiri- 

cally in the current version of the method. This formulation permits a 
dynamic overshoot of Cl above its maximum static value, but not of Cp. 

Also, note that aj)g is always less than since Aa^ = 0 when a - 0. 


5.4.2 UTRC g. A, B method (Refs. 27, 28) - A table-lookup correlation 
method, which has recently been synthesized into a more sti.',imUned format, 
was developed at United Technologies Research Center to determine C| and Cv. 
from three independent pariuiieters of the airfoil motion. The choice of the ' 
parameters was inspired by thin-air foil potential theory; they are inst.antan- 
eous incidence n(t), angular velocity parameter A = nc/2U^, and angular 
acceleration parameter B = ac /4Uc„‘'. It is interesting to note that these 
three parameters are sufficient to define to engineering accuracy the matclilug 
a(t) histories for the two cases shown in Klg. 17. The data base for the 
empirical correlations came from expeirlments on an NACA 0012 airfoil oscillat- 
ing in pitch ut M < 0.;i, including both sinusoidal ,md nonslmisoidal motion. 
None of the dat.a attained the extreme values that have been observed to 
accompany the fully developed vortex-shedding phenomenon of deep stall. 


-‘^‘■Ay3 MIT method (Ref. 29) - This method is basically an empirical repre- 
sentation of the fill Cl’S and moments due to the vortex— shedding plienomeiien for 
ramp clianges in angle of attack. The .actual angle of dynamic stall must be 


■ f gii ol tlio ladic.'il tenii is taken to be oppi'site to tlie sign 


of 
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, 1 +3“ baa normally been used. For 

ppeclfled separately; the value below static stall. Starting at 

4s < » < »DS> 'fc® '>“®“ ®'® . tlirTl^anv with time over a speoi- 

«'- ODS. “L ‘hi ^ f® «S“o Jerries that dapond on do/U„ 

fled tine Intetvai. from 1"®! Tf'^thlo la'^attalned before o - i%ax- then 
at the instant of dynamlu stall. U thl exponentially* with pro- 

Cl and Cjj remain constant the static stall values are attained, 

assigned time constants „ on the downstroke, when the 

uS^tanS atauc aeotlon charaoterletlca are reattalned. 

s.a.4 Looh heed method (Kefe. phSt-lag^tlme'censtants and 

empirical modeling or ayn^c 4cr4entB into a ficticious effective 

piLh-rate-dependent. etall-angU C„ from 

angle of attack. This linear combination of a- number of sepa- 

statlc airfoil characterlst c these elements are assumed to be analo- 

rate dynamic stall elements. S°m treated elsewhere in the literature, 

gous to the lag in circulation buildup on a pitching air- 

such as leading-edge jets, the lag walls, fluctuating pressure 

foil in potential flow. "^*thrvo«ax lift due to leading- 

propagation in i“®'>"''®'“,'’°"'^^L,o'^oloients are modeled directly from 
edge vortices on delta • ^j^^ating airfoils. In this sense, the 

SSVarmor:“d~rf “eedom than any of the other., and Information 
from many sources has been utilized. 

ff"thrpr2r-ruut« -ic! me ie the only method which dletingoiehes 
between pitching and plunging motion. 

5.4.5 Time - delay methods (Refs. ^ ^ jj g^parate universal dlmen- 

is that each dynamic ® At/c. regardless of the time history 

sionlees time ^irt 'ha parLnIter a - 2U..t/c. which U 

of the motion; t may be Indie ial aerodynamic response in the 

iraS-floi%rgS:!‘ m conetruction Of rbe " 

'ir.rgio^ofTti:ir:4i:^ ~ 

--^-1 r»* rtmnl 


t. = to + and Climax occurs a 2 _ taken trom 

(line segments I in the <^t,, Ci Continues this trend, but 

unsteady potential theory. ^moves rearward along the locus of the static 

the aerodynamic center of pr-sura moves rea^ard ai^^g ^ 

curve of X versus a. Therefore. Cm i interval, 

ing to the relation Cm « (0-25 - Kep^C^ curing 

,, m ^ r dprreases by an empirical exponential 

After lift stall at t = , calculated from the same relation 

law with respect to time, whereas C^i is calcuiatea . . x, ^ 

^ - A.-. rt t-nV-m-l 
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Fig. 29, Sketch of the time~de7ay 
method. 


for the. contribution of the vortex- 
shedding phemomenon to Cj^ nnd Cm. 

The Weatlando voralon (Ref, 4) 

Includes a dynamic Icadlng-edgc stall 
criterion in Tj that partially 
accounts for compressihillty effects 
and changeovers from traillng~edge to 
leading-edge stall. This method also 
seems to bo the only one to recognize 
the simple drag relationship, 

Cj) ^ Cj^ tan n, that obtains once the 
leading-edge suction is lo^t at 

5.4.6 ONERA method (Ref. 331- 
In contrast to the preceding methods 
that essentially curve-fit experimen- 
tal data with various algebraic or 
transcendental functions, this recent 
method utilizes a system of ordinary 
differential equations. If F repre- 
sents a dynamic force coefficient, 

®.8' > Cl or C{.j, and Fq represents 
the corresponding static coefficient, 
then F is split into a "linear'* part 
Fj and a "nonlinear" part Fj ; the 
system is written as 


F = Fi + Fj (4) 

f*! + a^F^ = + ^3^ (5) 

dAFn 

Fj + a^F2 + agFj == -as'^'Fo + ^6« (6) 

Here AFq represents the departure of the static coefficient from its linear 
trend, as illustrated for lift in Fig. 30, t is dimensionless time, and 
, 32, . . ., aj are coefficients to be determined empirically. 


For illustrative purposes, sinusoidal pitch oscillations, a = + a, sin mt , 

are considered. Below stall, that is, < «ss » - 0 and Eq. (5) with 

constant coefficients closely approximates the behavior of unsteady, linear, 
thin-airfoil theory. In the stall domain, however, ap varies as a function 
o a and hence with time. Each Up is evaluated by local linearization 
about Up. That is, data from a series of small-amplitude oscillations about 
various values of are used to develop empirical relations for a,,(a„). 

This requires a large quantity of small-.implitudc data, but ONERA has devel- 
oped a systematic way of generating the.se data very efficiently, once the 
peculiar requirements of the method ate recognized. 


Figure 31 derived from Ref. 33 illustrates .some predictions for moderate 
amplitude oscillations using the a^'s derived from data for which a =1° 
or le.ss. Comparable re.sults were obtained for Uj. ’ 


f _t 1" 1 h od s - All the methods described above pro- 
vide estimates of the effects of unsteady incidence changes, but each lias 
some shortcomings. Each method manages to reproduce reasonably well most of 
tlie data sets that were used in Its development, but almost no compari.sons 
have been made between any given method and indi pendent .sets of data. 


One notable exception has bacn 
described in Refs. 2 to 4. With the 
assistance of the originators and 
their colleagues I the present writer 
was able to compare the results of 
several of the aforementioned methods 
to an independent deep stall experi- 
ment* Four predicted quantities were 
compared; the phase angle > uu, at 
the onset of moment stall; the phase 
angle and maximum value of the normal 
force coefficient; and the value of 

^Mrain ‘ 

The Boeing-Vertol and UTRC 
methods tended to underpredict the 
effect of the dynamic stall vortex on 
maximum force and moment from the 
experiment, especially for large a^, 
whereas the MIT and Lockheed methods 
overemphasized the importance of 
vortex shedding on Cm- The Sikorsky 
time-delay method gave reasonable 
values for and However, 

this was somewhat fortuitous because it tended to predict lift stall later 
than the experiment, thereby compensating for the failure to include the 
extra lift due to the vortex. The phase-angle results for the other methods 

were mixed. 

A preliminary examination of the newer ONERA method has been encouraging, 
but it has also revealed that it is not well suited to data sets that were 
geLrL^d vithout regard to Its amall-smplltude and hlgh-Eroquency requirements. 



Fig. 30. Sketch of the oomponents 
of the ONERA method. 



Burprlaingly large-amplitude reaulta were obtained 
at low reduced frequenciea, but not when both and k were relativelv 

large, e.g., - 10 and k « 0.10. Also, it does not appear to bo able to 

model changea in atall behavior, auch aa, for example, static traillng-edL 
converting to dynamic leading-edge stall trailing edge 


6. CONCLUDING RKMARKS 


rhe general features of dynamic stall on oscillating air Tolls have been 
explained in terms of the vortex-shedding phenomenon, and the important differ* 

knowledge has been enhanced somewhat by 

lotential theory, boundary-layer theory, and Navier-Stokes calculations How- 
ever, most of what we know today has come from experiments and seralempirical 
modeling of the fluid dynamic forces and moments. The prediction techniques 
that have been developed for helicopter applications arc considerably better 

•, t-r"" ® satisfactory results await 

d better theoretical understanding of the effects of various types of airfoil 
motion and three-dimensional effects. 

The calculation of the flow around an oscillating airfoil at high inci- 
dence is a formidable task. The qualitative features of dynamic stall air- 
furtf captured by the empirical prediction methods, but all of them need 
further quantitative Improvements. At the present time, they are used for ' 
engineering purposes in preference to the analytical and numerical techniques 
1 this is likely to c.ange, especially for the 

i domain, as the current efforts to develop analytical tools for 

mature. Experiments, theoretical modeling, and 
H ^ can all be expected to play major roles in assessing and 

P c ing the three-dimensional elfects encountered in practice, but which 
are poorly nmlerstood at the presoiit time. 
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